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ABSTRACT 

We have applied the gaussian auxiliary field method introduced by Mazenko to a 
non-conserved scalar system with attractive long-range interactions. This study 
provides a test-bed for the approach and shows some of the difficulties encountered 
in constructing a closed theory for the pair correlation function. The equation 
obtained for the equal-time two-point correlation function is studied in the limiting 
cases of small and large values of the scaling variable. A Porod regime at short 
distance and an asymptotic power-law decay at large distance are obtained. The 
theory, is not, however, consistent with the expected growth-law, and attempts to 
retrieve the correct growth lead to inconsistencies. These results indicate a failure 
of the gaussian assumption (at least in the form in which we use it) for this system. 



I. INTRODUCTION 

The phase ordering dynamics of systems quenched from the disordered phase to the 
ordered phase has been extensively studied [Q]. There is a general consensus that 
at the late stages of domain coarsening these systems enter a scaling regime 0, in 
which the equal-time, two-point correlation function has the scaling form 

C{v, t) ^ (0(x, t)0(x + r, t)) = f{r/L{t)) , (1) 

where (f) is the scalar order-parameter field, L{t) is the characteristic length scale at 
time t after the the quench, / is a scaling function, and angled brackets indicate an 
average over initial conditions (and thermal noise, if present). 

A first principles calculation of the scaling function has proved to be a most 
difficult task. Even for the simplest model dynamics, that of a nonconserved order 
parameter (model A) with purely short-ranged (SR) interactions, exact results 
are rare and available only for cases of limited physical interest [^] . 

In the past few years closed-approximation schemes for the two-point correlation 
function of the SR model A (SRMA) have been proposed by a number of authors 
^-[|n|], based on a mapping 4>{r,t) = 0(m(r,t)) between the order-parameter and 
an auxiliary field m{r,t) which has, near a domain wall, the physical interpretation 
of a coordinate normal to the wall. With this new variable the problem of describ- 
ing the field at each instant of time is transformed into a problem of describing 
the evolution and statistics of the wall network. This approach enables the use of 
a physically plausible and mathematically convenient gaussian distribution for m. 
Such a distribution is unacceptable for the order parameter field itself, since this is 
effectively discontinuous at the domain size scale. 

The application of this sort of approach to both non-conserved and conserved 
(model B) dynamics, with purely SR interactions, has recently received a criti- 
cal review by Yeung et al. |]12[. Methods based on a description of the wall 
dynamics lead to an approximate Hnear equation for m(r, t), or for its correlator 
(m(x, t)m(x + r, t)), A different and promising approach, due to Mazenko [0, 
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aims at deriving a closed n on- linear equation for C(r, t), built on the equation of mo- 
tion for model A, using the single assumption that the field m is gaussian distributed 
at all times. It has the advantage of yielding results with a non-trivial dependence 
on the spatial dimension d and is also easily extensible to 0{n) component systems 
with topological defects, i.e. with n < d |10[ . 

The only uncontrolled feature of this approach is the gaussian assumption. Re- 
cent simulation tests have shown, however, that this assumption is not entirely 
satisfactory: Blundell et al. |]13| have made an absolute test (free of adjustable pa- 
rameters) of the relation between two different scaling functions, revealing a disap- 
pointing agreement with the theory. The discrepancy decreases, however, in higher 
dimensions, in agreement with an argument that the gaussian approximation 
becomes exact in the limit d oo. Yeung et al. [Q, using data of Shinozaki and 



Oono IjT^ for = 3, have checked the single-point probability distribution for m, 
finding it to be flatter at the origin than a gaussian. It is not difficult to derive an 
analytical expression for the two-point distribution P(m(l), m(2)), valid for 
m(2) and |r| small compared to L{t) |T5|. It differs from a gaussian for fixed spatial 
dimension d, but is consistent with a gaussian in the limit d ^ oo. 

Despite these reservations, the gaussian approach has been shown to give good 
results for the SRMA, displaying most of the expected physical properties 0, |TU|. 
In this paper we try to extend the limits of this approach by applying it to model A 
dynamics with attractive long-range (LR) interactions. This application addresses 
a basic difficulty, not necessarily caused by the use of the gaussian assumption: 
the attempts to extend the approach beyond the simplicity of the SRMA produce 
equations for C(r, t) which do not seem to respect the expected growth- law for the 
typical domain size L{t). A lack of a proper scaling of the terms in the equation for 
C(r, t), derived naively, is apparent for the case of a scalar order parameter, namely 
for the LR model A (LRMA) and for the SR model B (SRMB), although not for a 
vector order-parameter in which case a 'naive' dimensional analysis of the equation 
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agrees with the known growth law (with logarithmic corrections for n = 2) [|T^. We 
shall see how this situation arises for the LRMA and present our understanding of 
it. In the case of the SRMB, a naive application of the method, however, omits the 
important bulk diffusion process which plays a vital role in the coarsening. Mazenko 
10] has attempted to solve the problem by accounting explicitly for the bulk diffusion. 
It is not clear, however, that any analogous mechanism is present here. 

The scalar case, which is usually the more interesting one in the applications, is 
exceptional because an extra length, time independent at late stages, the domain- 
wall thickness, plays a role in the dynamics, and therefore power counting of lengths 
by dimensional analysis may not yield the right scaling (in terms of the characteristic 



length) for the different parts in the equation of motion. For the SRMB [|T6l and 
the LRMA fl^, |l8l dynamics the growth laws are L{t) ~ t^/^ and L{t) ~ ^1/(1+'^) 
(for cr<l), respectively, for n = 1, and L{t) ~ t^/^ and L{t) ~ respectively, for 



n > 2 (with logarithmic corrections for n = 2 fl^); where 0<cr<2 is the exponent 
describing the LR interactions, which decay as . For n = 1 and l<a<2, the 

long-range interactions are irrelevant and the growth-law is the same as for the SR 
case |TS[. The SRMA, however, is exceptional since the predicted growth-law. 
Lit) ~ t^/^, is the same for both the scalar and vector order-parameters, accidently 
allowing for a 'naive' dimensional analysis of the scalar equation of motion to agree 
with the growth-law. In this case the role of the extra length in the scalar equation 



can be ignored as the result of two canceling errors |T8[. Therefore we wonder 
if the success of the Mazenko method with this scalar model might be somewhat 
fortuitous. In other words, we raise the question of whether this approach (or 
any other closed theory), naively applied, can succeed for those dynamical models 
where naive dimensional analysis gives the wrong growth law. In this respect it is 
interesting that a straightforward application of the method of Kawasaki, Yalabik 
and Gunton to the LRMA also gives the wrong growth law for n = 1, i.e. it 
gives the t^/'^ growth suggested by naive dimensional analysis. 
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In this paper we have developed an extension of Mazenko's approach for the 
LRMA. Besides having interest by its own right, the study of this model provides a 
test-bed for the approach and shows some of the difficulties any approximate closed 
theory must resolve. 

II. THE MODEL WITH LONG-RANGE INTERACTIONS 

We consider a system with long-ranged attractive interactions, falling off with dis- 
tance as r^(d.+(T) ^ suitable Hamiltonian functional of the scalar field is 

H[<p]= Jd'r[{V<pf/2 + V{<P)] + {JLR/2) Jd'^r |rfV[0(r) - 0(r')]V|r - (2) 

where as usual we have taken the short-range part to have the Ginzburg-Landau 
form, Ji/j > 0, and V{(f)) has a local maximum at = and global minima at 
= ±1. The model is well defined for < cr < 2. The equation of motion for a 
non-conserved field reads dip/dt = —SH/6(j), i.e. 

^^ = V'<P-V'{<j>) + VU<P), (3) 
where V'{(j)) = dV/dcj) and the LR force is given, both in real and Fourier space, as 

ViM) = JLRjd'r'W)-<P{v)]l\r-v'r^ (4) 

d'^k 



and 



Jlu h{d, a) I ^ 0(k) re^^-'^ , (5) 



hid,a) = Qid,^)^^!^y (6) 



In noise is absent since temperature is an irrelevant variable [20|. From an 
analysis of (|^), assuming the validity of the scaling hypothesis (1), the following 
growth-law has been predicted for a scalar order-parameter |18| 

L{t) ~ , 0<a<l 

~ , l<a<2 , (8) 



in which the crossover a = 1 separates the regime where domain growth is faster due 



to the LR correlations from the regime where these become irrelevant [18], ^1 



III. THE SCALING EQUATION 

To obtain an equation for the two-point correlation function (|1]) we multiply (^), 
evaluated at point (1) = (ri,ti), by evaluated at point (2) = (r2,t2) and average 
over the ensemble of initial conditions yielding, at equal-times, 

iac(i,2) 



2 dt 



V'C{1, 2) - (0(2)\/'(0(l))) + mvLnim)) ■ (9) 



We will call (0(2)\/'(0(l))) and (0(2)V^^(0(1))) the 'non-linear' (NL) and the 'long- 
range' (LR) terms of the equation for C(r, t), where r = ri — r2. In (P) the LR term 
reads, both in real and Fourier space, 

mvUHl))) = JLRjd'r'[C{r',t)-C{r,t)]/\r-rr^'^ (10) 
= Jlr h{d, <y)j^, C^(k, t) k'^e^'^ , (11) 

Assuming the existence of a late-time single-scaling regime, we expect C(r, t) to 
take the scaling form (|l|), in terms of which (|^) reads 

- \^ = j^if" + - {mv'm))) + {mvu<f>m , (12) 

where x = r/L{t) is the scaling variable and /' = df/dx, etc. In the equation above 
L/L 1/t, if L{t) grows as a power-law. The LR term now reads 



imvum))) = -j;i^jd'x'[fix')~f{x)]/\^-^r^- (13) 

JLRh{d,a) f d'^y 



/t^ . (14) 



L{tY J (27r)^ 

where g{y) is the Fourier transform of /(x), and y = kL{t). 

From an analysis of (^) for an isolated, stationary, planar wall, we find that to 
leading order the equilibrium planar wall profile saturates as 



l-02(r)~-^, (r-^oo), (15) 



where Vq = {(PV/d(j)'^)^2^i and r is the distance from the walL Hence we expect that 
throughout the bulk region |0| will be below saturation by an amount ~ l/L{tY . 
Even with this power-law decay we still expect there to be well-defined walls, with 
a time-independent 'thickness' defined for example from (15) via w" = JLi?,/K)". 
Therefore, domain walls may be regarded as 'sharp' at late-times, when L(t) ^ w. 
It follows that Porod's law [Q, g{y) ~ A{d,cr)/y'^^^ for y ^ 1, holds within the 
regime kw <^ 1 <^ kL(t) = y (corresponding to w <^ r <^ L(t) in real space), 
in which case eq. (0) yields for the leading scaling behaviour of the LR term, as 
X ^ 0, 



J LRh{d,a) / r d'^y ^ A{d,a) h{ d,l-a 

L{tY [J {27iY^^^^^ (27r 

= ^^^(^^f'"^ f^^^'")^^f;'""^^ + 0(l)V l<a<2. 

This result will be exploited below to determine the amplitude A{d, a) of the Porod 
tail, within the gaussian approximation. 

IV. THE GAUSSIAN APPROXIMATION 

In order to transform (^ or (|T2]) into a closed equation we need to express the 
NL term as some approximate non-linear function of C{r,t). A key idea, exploited 
by several authors [6-11] within SR model A dynamics, is to employ a non-linear 
mapping between the order parameter 0(r, t), which at the scale of L{t) is effectively 
discontinuous near walls, and an auxiliary 'smooth' field m{r,t), whose zeros define 
the wall network. This introduces the wall structure into the problem and allows 
the approximation to be implemented through the new field. 

From the equation of motion (^ we can see that, just like in the SRMA, if the 
initial field satisfies |0| < 1 then this condition will hold at all times, assuring that 
a one-to-one mapping can be defined. For this model we have in mind, following 
Mazenko's treatment for SR interactions, to identify the field m(r, t) at points r near 
domain walls as the (signed) distance to the nearest wall (along its local normal), 
with the sign of m being that of 0. This determines m uniquely when m <^ L{t). 
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To specify m everywhere in space, we define tlie function <^(r, t) = (j){m{r,t)) by 
extending Mazenko's suggestion [0 of using the equihbrium planar domain wall 
profile function for an isolated wall, with m the coordinate normal to the wall, i.e. 
the function 0(m) is specified by the equation 







[{m' — niY + y^] 2 



(17) 



with boundary conditions 0(0) = and 0(m) — sign (m) for |m| oo. Using (|T 
we rewrite the NL term in (|9|)- ([l2|) as 

rfV(^(l)) 



{mv'm))) = {^m- 

Q{d,a) 



dm{iy 



+ 



+ 00 



ds' 



LR 



ll+cr 



(0(m(2))[0(m(l) + s)-0(m(l))]) (1^ 

im(l)k+im{2)k' \ 



where h{d, a) and Q{d, a) are given by (P)-(0), and (f){k) is the Fourier transform of 
(f){m). 

Following Mazenko 0], we now make the key assumption that m(r, t) is a gaus- 
sian field (with zero mean) at all times, with a pair distribution function 

1 



P(m(l),m(2)) = N exp 



^o(l) 



m(l)2 m(2)2 27m(l)m(2) 
2(1-7^) ^:VTy+ 5^27" ^5o(l)5o(2)^ 

1 



m(l)n, ^(1,2) = iV 



(20) 



(21) 



/So{l)So{2) 27r^(l-72)^o(l)5o(2) 
We also note that, as the walls become effectively sharp in the late-time regime, 
we can use the profile 0(m) = sign (m) to evaluate the leading contribution to the 
scaling functions. From ( plSf ) we expect the effect of ignoring the power tail in the 
profile is to neglect a quantity of relative order ~ 1/L{ty in the LR part of (0). 
The purely SR part of the NL term is then simply given, as a non-linear function of 
C{r,t), by Mazenko's result for the SRMA 

d^(f){m{l))\ 2 



0(2)- 



dm{l] 



(22) 
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Deriving a similar result for the LR part of the NL term is more tricky. There 
are three different ways to perform the calculation: we will outline the basic steps 
of each one. Representing 0(m) in Fourier space and Taylor expanding 0(m + s) 
in powers of s, using the gaussian property and returning to real space, gives the 
formal expansion |23| 



/+00 / 
dsy^ (0(m(2))[0(m(l) + s) - 0(m(l) 

= Q{d, a) dsj^ E ^2 . (23) 

Using C(l,2) = (sign (m(l))sign (m(2))), the integral representation sign (m) = 
l/(z7r) J^^dzexp[izm]/z, and the gaussian property, the series can be summed. 



Finally, differentiating with respect to Co(l,2) = 5'o(l)5'o(2), performing the z 
and s-integrals and integrating back, yields the non-linear function 

FNL{h2) = -^^liyn //''^^^^^^W ' (24) 

a{d, a) = h{d, a) ^7^^ . (25) 

Alternatively, we can take 0(m) = sign (m) from the start and do the s-integral, 
giving 

P ^^ 2JLRQid, a) / sign (m(l)) sign (m(2)) \ 

F^L 1,2 = ( r^TTT^ /' 26 

use integral representations for signm and l/lml'^, do the gaussian integral, differen- 
tiate with respect to Co(l, 2), perform the remaining integrals, and finally integrate 
back over Co (12), yielding the same result. Taking into account (|l^) and using a 
signm profile, (^) can be recognized as the leading order result for a 0^-potential 
NL term, i.e. < 0(2)0(1) (1 — 0^(1)) > Finally, the simplest derivation is to take 
the integral representation of signm in (p6D, use the gaussian property, differentiate 
with respect to Co(l,2) and do the gaussian integral, leading to the same point as 
the first calculation before its final integrations. This derivation, however, does not 
provide the appealing intermediate expressions (|23|) and (E6|). 
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According to our identification of m{r,t) as a distance from tfie interface, we 
expect 5*0 = (m^) to have the scahng form So = L{t)^, which can be used along with 



I) and (|24D to rewrite equation (0) for the scahng function in the form 



I L ^, 1 / d — 1 ^, 2 / TT ^ 

^ ■ '''' 

For a < 2 the SR part in (P?]), scahng as is neghgible compared to the LR 

part, scahng as l/L", and can be ignored (but see the discussion in section VI!). 
Demanding that the left side of (^T]) balance the terms of order on the right 
requires L/L '^/L'^i i-e. that L{t) ~ Note that this disagrees with the 

expected form (||)! In section V we will argue that a resolution of this discrepancy 
requires us to drop the left side of p7| ) in leading order. For the moment, however, 
we pursue the original (and a priori natural) assumption that the left side scales as 
and write 

m = {jLRi^tf" , (28) 

where /i is to be determined. Dropping the SR terms from (|27|) gives the final 
equation for the scaling function /(x): 

= {fi/2a) xf + f " ^) ^^^"^^) • (29) 

Equation (^) has to be solved numerically for general scaling variable x. How- 
ever, it is straightforward to derive analytically the behaviour for small and large 
X. Using the Porod's law form f{x) = 1 — a{d, cr)x + ... for small x (it is simple to 
show that this is only consistent short-distance behavior), we find that the LR part 
of the NL term, (24), has a leading scaling behaviour as x which is similar to 



m 



where B{x,y) is the beta function. Performing a small-x expansion of eq. (p9|), we 
find that the dominant terms for x — are obtained from the terms multiplying 
Jlh in (27), whose small-x expansions are given by (p!6D and (|30D. Matching powers 
of X for general 0<cr<2, and using A{d,a) = —a{d,a){2TTY/h{d, 1) in ([16]) (which 
follows from Fourier transforming the Porod tail [^), we find 

«(^>^) = V [ri&y) ' o<^<2' (31) 

for the coefficient of x in the small-x expansion of f{x) |2^. For a = 2, this reduces 



to the SR resuh a{d, 2) = 2/{n^d- 1). 

For 0" < 1, ( PH ) was obtained by matching the terms of 0{x^^'^) in (p!6| ) and (pO|). 
The leading (constant terms) yield an interesting sum rule to be satisfied by the 
structure factor scaling function g{y) for this range of a: 
d'^v 2'^/^ 

— ^ g(y) y'' = < a < 1 . (32) 

{2-kY r(^) sin(i±^7r) ' ^ ' 



We now look at the large-x asymptotic form of eq. (p9D , and discuss the large-x 



behaviour of f{x). In this limit, f{x) and the final two terms in (^) become 
g{Q)/x'^^'^ and —a{d, a) Txf{x)/2 respectively. In this regime, (29) can be integrated 
to give 

fix) ^^^(0) ±- + -J— (33) 

where we note from (25) and (6) that a is negative. In general, both terms in (33) will 
be present in the large-x solution. On physical grounds, however, we do not expect 
f{x) to fall off with distance more slowly (in a power-law sense) than the underlying 
interactions, which decay as r~(^+'^^ (An exception is when sufficiently long-range 
power-law spatial correlations are present in the t = state. This power-law can 
then persist for general times Here, however, we consider only short-range 

correlations in the initial state.) We infer that either /x < \a\TTa / {d + cr), so that 
the second term in (33) is subdominant for large x, or A = 0. The first possibility, 
however, implies that the coefficient of the (dominant) first term in (33) is negative 
(since g{Q) > by definition), i.e. f{x) approaches zero from below, which also seems 
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unphysical (and disagrees with numerical simulations [0). We conclude that the 
only physically sensible possibility is that A vanishes in (33). This can, presumably, 
only happen for a special choice of /i, so the condition A = determines fi. This 
mechanism is very similar to that which determines fi for short-range interactions 
0, |10[- Note that, if /(x) is to approach zero from above for x ^ oo, (33) gives the 
inequality 

/i > 7r\a\/{d + a) . (34) 

A sum rule for fi can be obtained by integrating (29) over space: 

2a\a\ f , rf/"^ 

^ = —±i / d'^x / de sec" e . 



dg{0) J Jo 

Finally, it should be noted that the above analysis implicitly assumes that g{0) 
is finite, i.e. that f{x) decays faster than x~'^. In fact, the mathematical structure 
allows for f{x) ~ x'^ with p < d [^, but we reject this possibility on the physical 
grounds that we appealed to before, namely that, at least for initial states with 
only short-ranged spatial correlations, the scaling function should not decay with a 
smaller power than the underlying interactions. 

V. TWO-TIME CORRELATIONS 

The gaussian approach can also be used to evaluate the two-time correlation function 
C(r, ti,t2) = (0(x, ti)0(x + r, ^2)) and, in particular, the autocorrelation function 
^(^1)^2) = C{0,ti,t2). The calculation is simplest in the limit ^2 ^ ti, when C ^ 
and the full nonlinear equation can be linearized, Fourier transformed, and explicitly 
integrated. In this regime the analog of (9) for two-time correlations reads, in Fourier 
space (dropping the SR term on the right), 

^ = -JLR\hrC^ + ^C^, (35) 

Ot2 2 fit 

where (28) has been used for L{t). We integrate (35) forward from time ati, where 
a 3> 1 ensures that the condition ^2 ^ ^1, required for the validity of (35), holds at 
all times. This gives 

a(ti, t2) = Ck(ti, ah) exp{- JiR|/i|r(t2 - ati)} . (36) 



at 
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Using the scaling form Ci^{ti,ati) = Lfga{kLi), where Li = L{ti), and summing 
over k for t2 ^ ati gives the autocorrelation function 

A{ti,t2) = const (Li/Ls)''-"''^!'^/''^ , (37) 

where 'const' is clearly independent of a. The physical requirement that A decrease 
with increasing ^2 gives the inequality fi > 7i\a\a /2d, which is guaranteed by (34) 
for d> a. 

The connection between the parameter /i and the exponent describing the decay 
of the autocorrelation function is similar to that obtained for purely short-ranged 
interactions [[l^, ^ . 



VI. DISCUSSION AND SUMMARY 

We have extended the original Mazenko gaussian approach to the LR model and 
evaluated the late-time leading contribution to the NL term of eq. (0), yielding a 
dominant LR part given by (p^, which is of order 1/L" . An infinitely sharp wall 
profile has been used which amounts to neglecting a quantity of relative order 
The LR term in the equation, (|11])-(|T^, is of the same order and has an amplitude 



which is a function of x, d and a, but its non-local nature (i.e. its dependence on 
the values of f{x) everywhere) makes the problem particularly hard to handle. 

Despite the profile power-law decay ([T5|) induced by the LR interactions, the 
scaling function exhibits Porod's law, i.e. a linear short- distance behaviour in real 
space with coefficient given by (0). This is consistent with the assumption that 
at late-times there are well-defined walls with a constant 'width' independent of 
L{t). This is an important point of principle, on which the identification of the field 
m(r, t) and the mapping (|T^ rely, and also a key ingredient in the first-principles 



derivation |T^, |18[ of the growth-law (j^). 

The central question we want to address in this paper is whether the gaussian 
theory is able to yield the correct growth-law for this model. We have seen that the 
'naive' application of the gaussian approach present in section IV ostensibly gives 



the wrong growth law: (pSf) instead of (§). A related problem is the SRMB to which 
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Mazenko has attempted to apply the gaussian approach , yet the correct growth- 
law does not come out of the theory as cleanly as in the SRMA. In this system local 
conservation imposes a bulk diffusion process which controls the interface motion 
and delays domain coarsening relative to the purely relaxational dynamics of model 
A. There are some common features between the dynamics of a conserved and a LR 
interacting field, namely the existence of a bulk profile which relaxes rapidly to a 
non-saturating value as the walls move. One key difference, though, is that the true 
growth law for SRMB (t^^^) is faster than that obtained by a naive application of 
the gaussian approach (without allowing for bulk diffusion), which gives t^^^. 

Before implementing any approximation we focus the analysis on the exact equa- 
tion (p!2D . If the growth-law (H) holds, the time-derivative term must be negligible 
compared to the LR term ([T^) , which scales as L~'^, and therefore the NL term must 
have a leading contribution of order which exactly cancels the LR term in the 
scaling limit. In fact, this condition determines the late-time leading contribution to 
the scaling function. Within the gaussian approximation, it amounts to neglecting 
the first term in ( [29| ) (which came from the left-hand side of (0)), to give 



Solving this equation gives the scaling function /(x), within the gaussian approach, 
provided the growth law is slower that t^^'^. However, there seems to be no way 
to determine the growth law within this scheme. Moreover, ( |38D has a serious 
shortcoming. If we integrate the equation over all space, the first term drops out, 
giving the sum rule 



Since the integrand is positive definite, the only way this sum rule can be satisfied 
is for f{x) to be negative for some range (or ranges) of x, with sufficient negative 




(38) 




(39) 



weight to satisfy (P5|). This seems a priori improbable for a nonconserved order 



parameter, and indeed numerical simulations ||2^ show no hint of it. 
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We emphasize, however, that since our fundamental equation (|T^) is exact, the 
analogue of (36) obtained without making the gaussian approximation must be ex- 
actly true. Because the true growth is slower than t^^" , the left side of (12) is 
negligible in the scaling limit. Taking the Fourier transform of the equation, and 
setting A; = 0, the 'long-range' term vanishes. This leaves the identity 

d'x (0(2)l^'(0(l))) = , (40) 

of which (39) is the special case obtained within the gaussian approximation. 

Our results seem to indicate that the gaussian approach, applied to the bulk 
equation of motion, is unable to account for the qualitative feature of coarsening in 
systems with long-range interactions. However, we cannot make a definitive state- 
ment as we have not exhausted all the possible choices for the gaussian field. There 
may exist a mapping definition which is physically more appropriate and works 
better then (p!7|). What seem to be clear is that, beyond the simple nonconserved 
system with short-range interactions, one cannot apply the gaussian approach in 
a straightforward and 'naive' manner to construct a closed equation for the scal- 
ing function. Just as for the conserved scalar system, a deeper understanding of the 
underlying physics may be required in order to implement a more controlled approx- 
imate scheme. It is possible that this might be achieved by means of an interfacial 
approach [^ . 

To summarize, the gaussian approach, naively applied, is unable to yield a 
growth-law different from that obtained from dimensional analysis of the linear 
terms in the equation of motion. As mentioned above, the failure of the gaussian 
approach in this context could be due to the particular choice employed for the 
mapping between and the gaussian field. For example, one can in principle use 
the same mapping as Mazenko [0, 0"(m) = K'(0(m)), appropriate to a purely 
short-range interactions. However, this leads to an inconsistent scaling analysis of 
the equation for C(r,t) (e.g. at short-distances there is no LR part in the NL term 
to match the LR term (p!6|) ), and m(r,t) can no longer be regarded as a distance 
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from a wall. By contrast, the mapping employed here, defined by eq. (|T7D, seems 
far more natural and physically suitable for a system with LR interactions. 

Finally we note that the present methods can also be used for a vector order pa- 
rameter with long-range interactions. In that case the t^^" growth obtained within 
the gaussian approach is correct (apart from logarithmic corrections for n = 2 
The purpose of the present paper, however, is to test the method on those systems 
which provide the greatest challenge, i.e. scalar systems, in the hope that the diffi- 
culties identified here may stimulate the development of more robust approximation 
schemes. 
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